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THE  DISTRIBUTIONS  OF  MEASURES  OF  VISIBILITY 
ON  LINE  SEGMENTS  IN  THREE  DIMENSIONAL 
SPACES  UNDER  POISSON  SHADOWING  PROCESSES 


by 


M.  Yadin  and  S.  Zacks 


Abstract 


The  problem  of  determining  the  moments  of  visibiltiy  measures  on  star-shaped 
curves,  under  Poisson  random  shadowing  process  in  the  plane  (see  [ 4  ])  is 
generalized  here  to  the  problem  of  determining  the  moments  of  visibility 
measures  on  line  segments  in  a  three-dimensional  space,  when  the  shadowing 
process  is  created  by  a  Poisson  random  field  of  spheres  of  random  size.  The 
present  paper  provides  the  methodology  for  reducing  the  three-dimensional 
problem  to  a  two-dimensional  one,  and  employing  the  general  theory  of  [4  ] 
to  solve  the  present  problem.  Although  the  general  approach  is  similar,  the 
functions  derived  in  the  present  paper  are  somewhat  different  for  the  sake 
of  simplifying  the  computations. 

Key  Words;  Poisson  Shadowing  process;  Lines  oh  Sight;  Visibility  Pfiobabltltles 
Measunes  oh  Visibility;  Moment i  oh  Visibility. 
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1.  Introduction 

Consider  a  source  of  light  located  in  the  origin  and  an  arbitrary 

3 

straight  line  segment)  C  ,  in  a  three  dimensional  Euclidean  space,  R 
A  layer  of  random  spheres  is  located  between  the  origin  and  C  .  The 
centers  of  the  spheres  are  randomly  located  between  two  parallel  planes, 
according  to  a  Poisson  process  with  a  specified  intensity.  Furthermore, 
the  radii  of  the  spheres  are  independent  identically  distributed  random 
variables,  having  a  specified  distribution,  such  that  both  the  origin, 

0  ,  and  the  segment,  C  ,  are  not  covered  by  the  random  spheres.  A 

point  P  on  C  is  called  visible  (in  the  light),  if  the  line  segment 
OP  does  not  intersect  any  random  sphere.  The  measure  of  visibility  on 
C  is  the  total  portion  of  C  which  is  visible.  We  are  interested  in 
the  distribution  of  this  random  variable. 

In  a  previous  study  [ 4  ]  we  developed  the  methodology  of  determining 
the  moments  of  the  measure  of  visibility  on  star-shaped  curves  in  the 
plane  and  provided  also  an  approximation  to  its  distribution.  The  objec¬ 
tive  of  the  present  paper  is  to  develop  a  similar  methodology  for  treating 
three  dimensional  problems.  For  this  purpose  we  devote  Section  2  to  the 
discussion  of  the  geometry  of  the  three  dimensional  shadowing  model.  The 
three  dimensional  problem  is  reduced  in  Section  3  to  a  two  dimensional 
problem.  We  distinguish  between  two  cases:  the  case  where  the  plane  con¬ 
taining  0  and  C  is  perpendicular  to  the  layer  of  the  random  spheres 
and  the  case  where  that  plane  is  slanted  with  respect  to  that  layer.  In 
the  first  case,  a  standard  Poisson  field  of  spheres  (uniformly  distributed 
centers  with  radii  independent  of  the  location  of  the  centers)  is  reduced 
to  a  standard  Poisson  field  of  (intersecting)  disks  on  the  plane.  This  is 
not  the  case,  however,  when  the  plane  is  slanted.  Subsections  3.1  and  3.2 
are  devoted  to  these  two  cases,  the  standard  and  the  non-standard  ones.  In 
Section  4  we  introduce  the  notion  of  visibility  measure  and  provide  recur¬ 
sive  formulae  for  computing  its  moments.  Section  5  provides  a  development 
of  certain  basic  functions,  called  here  the  K-f unctions,  required  for  the 
determination  of  visibility  probabilities  and  for  the  moments  of  visibility 
measure.  In  Section  6  we  further  derive  the  K-f unctions  under  the  assump¬ 
tion  that  the  radii  of  the  spheres  are  uniformly  distributed  over  the  inter¬ 
val  [0,b]  .  We  probide  also  some  numerical  illustrations  of  the  standardized 
moments  of  the  measure  of  visibility.  Derivation  of  the  first  two  moments  of 


the  distribution  of  the  radii  of  disks  in  the  plane  intersecting  the 
layer  of  spheres,  and  derivation  of  the  corresponding  non-central 
moments  is  provided  in  appendices.  The  theory  provided  here  has  a 
variety  of  important  applications.  For  example,  a  helicopter  is 
flying  over  a  forest,  and  its  flight  path,  at  a  specified  time  period, 
can  be  described  by  a  given  vector.  An  observer  is  located  at  a  given 
point  in  the  forest.  The  crowns  of  trees  are  located  at  different 
heights,  and  can  be  modeled  (approximately)  as  spheres  of  random  radii. 
The  problem  is  to  determine  the  distribution  of  the  visible  portion  of 
the  flight  path.  Instead  of  helicopter  and  trees  in  the  forest,  one  can 
think  of  an  airplane  and  a  layer  of  clouds.  There  are  many  other  appli¬ 
cations  in  areas  of  communication  with  sattelites,  astronomy,  biology, 
fishery,  etc. 


2 .  The  Three  Dimensional  Shadowing  Model 


Let  U  ,  W  and  C  be  three  parallel  planes,  having  distances 

it  it  it 

u  ,  w  and  r  ,  respectively,  from  a  point  0  .  It  is  assumed  that 

0  <  b  £  u*  <  w  Sr  -b  ,  (2.1) 

where  b  is  some  specified  positive  value.  Let  C  be  a  straight  line 
in  C  and  A1  the  plane  passing  through  0  and  C  .  Let  W  and  U 

ft  fg 

be  the  straight  line  in  which  M  intersects  W  and  U  ,  respectively. 

The  three  parallel  lines  U  ,  W  and  C  are  at  distances  u,  w  and  r 
from  0  ,  satisfying  the  relations 


* 

r 

—  «*  cos  d> 
r 


(2.2) 


it 

where  «p  is  the  angle  between  the  axes  Z  and  Z  (see  Fig.  1)  . 

it  ft 

Notice  that  Z  is  perpendicular  to  C  ,  while  Z  lies  on  M  and  is 

* 

perpendicular  to  C  and  both  Z  and  Z  pass  through  0  . 

A  source  of  light  is  located  at  0  and  spheres  of  random  diameters 

are  randomly  centered  between  U  and  W  .  It  is  assumed  that  the  centers 

of  the  spheres  follow  a  homogeneous  Poisson  process  of  intensity  X  , 

*  * 

between  U  and  W  .  Furthermore,  the  radii  of  the  spheres  are  independent, 

identically  distributed  (iid)  random  variables,  X,,,  ...  having  a  common 

c.d.f .  F(x)  concentrated  on  [0,b]  .  Condition  (2.1)  implies  that  none  of 

these  random  spheres  either  covers  0  or  intersects  C*  .  These  spheres 
* 

cast  shadows  on  C  and  in  particular  on  C  . 

For  the  purpose  of  studying  the  visibility  on  the  line  C  we  consider  the 
two-dimensional  shadowing  process  on  M  .  Shadows  on  C  are  cast  by  random 
disks  on  M  ,  which  are  generated  by  random  spheres  (intersecting  M  ) . 


3.  Reduction  to  a  Two-Dimensional  Shadowing  Model 


In  the  present  section  we  consider  the  reduction  of  the  shadowing  model 

to  M  .  In  particular,  we  have  to  determine  the  distribution  of  the  centers 

of  shadowing  disks  on  M  and  the  conditional  distribution  of  their  radii 

given  the  location  of  their  centers.  In  Fig.  2  we  present  the  geometry  of 

the  reduction  to  a  two-dimensional  model  by  the  projection  of  the  three- 

dimensional  space  on  a  plane  perpendicular  to  C  .  We  restrict  attention 

to  spheres  which  are  centered  within  a  prism,  P  ,  whose  projection  is  the 

parallelogram  ABCD  (see  Fig.  2).  The  distances  from  M  of  spheres  centered 

outside  this  prism  is  greater  than  b  .  Accordingly,  such  spheres  do  not 

intersect  M  and  therefore  they  cannot  cast  shadows  on  C  .  A  sphere  centered 

within  the  prism  P  intersects  M  if  its  radius,  x  ,  is  greater  than  its 

distance  from  M  ,  d  .  The  intersection  of  such  a  sphere  with  M  is  a  disk 
2  2  1/2 

of  radius  y  -  (x  -d  )  .  For  simplicity,  we  will  assume  that  all  the 

spheres  centered  within  the  prism  P  generate  on  M  disks;  those  spheres 
with  x£d  generate  on  M  a  disk  with  radius  y«0  .  The  center  of  a  disk 
generated  by  a  sphere  is  the  projection  of  its  center  on  M  .  Such  projec¬ 
tions  may  lie  between  two  parallel  lines  on  M  ,  having  distances  u-&  and 
w+B  from  0  ,  where 

S  -  b  tan  *  .  (3.1) 

3.1  Random  Disks  on  M  When  8  ■  0 

In  the  case  of  3-0  ($-0)  the  plane  M  is  perpendicular  to. the  planes  C*  , 

(Jj  and  U  .  In  this  case  ^*  "  ^*  *  ~*  “  1  •  The  distribution  of  the 

r  w  u 

location  of  the  center  of  disks  on  M  follows  in  the  present  case  a  homo¬ 
geneous  Poisson  process,  with  intensity  y»2Ab  ,  between  two  parallel  lines 
U  and  W  on  M  ,  having  distances  u  and  w  from  0  .  Indeed,  consider 

any  Borel  set,  B  ,  on  M  ,  of  area  A  located  between  U  and  W  .  Disks 

which  are  centered  in  B  are  generated  by  spheres  centered  in  the  cylinder 
set  of  B  within  the  prism  P  .  The  volume  of  this  cylinder  set  is  2bA  . 
Accordingly,  the  number  N{B}  of  disks  on  M  centered  in  B  has  a  Poisson 
distribution  with  mean  2AbA  . 

Furthermore,  since  centers  of  spheres  in  P  are  uniformly  distributed,  the 
distances  of  these  centers  from  M  are  i.i.d.  random  variables,  having  a 
common  uniform  distribution  on  [0,b]  .  Accordingly,  the  c.d.f.  of  the 
radius  of  a  disk  on  M  ,  Y  ,  is 


(3.2) 


G(y)  =  £  /  F((y2+u2)1/2)  du 


In  the  special  case  that  X  has  a  uniform  distribution  on  [0,b]  we 


obtain 


G(y)  -  1  -  \  (1-(£)V/2  -  |  (£)2  log  (J) 
+  |  (£)2  log  (1  +  (1-(£)2)1/2) 


0£y£b 


(3.3) 


Notice  that  in  this  uniform  case  G(0)  ■  1/2  .  It  is  expected  therefore 
that  half  of  the  spheres  in  P  will  not  intersect  M  .  The  case  where 
the  Poisson  process  of  the  disks  centeres  is  homogeneous,  and  the  distri¬ 
bution  of  Y  is  independent  of  the  location  of  the  disks,  was  called  in 
[  4  ]  the  standard  case. 

3.2  Random  Disks  on  M  For  0  >  0 

Ti;e  projection  of  centers  of  random  spheres  on  M  ,  in  the  case  of 
0  >  0  ,  yields  centers  of  disks  which  follow  a  non-homo geneous  Poisson 
process.  Indeed,  the  volume  of  cylinder  sets  in  the  prism  P  ,  corres¬ 
ponding  to  Borel  sets  of  area  A  between  the  two  parallel  boundary  lines 
on  M  (of  distances  u-0  and  w+0  from  0)  is  not  constant.  We  can, 
however,  partition  the  prism  P  to  a  collection  of  thin  layers,  so  that 
the  centers  of  disks  generated  by  spheres  centered  in  each  one  of  those 
layers  practically  follow  a  homogeneous  Poisson  process.  More  specifically, 

ft  fc 

consider  a  layer  bounded  between  two  parallel  planes  V  and  \l  ,  having 

ft  ft  ft 

distances  v  and  v  +  Av  from  0  ,  where 


*  *  *  *  * 
u  S  v  <  v  +Av  <  w 

(see  Fig.  2).  These  planes  intersect  M  in  two  parallel  lines,  with 

* 

distances  from  0  of  v  and  v+Av  ,  respectively.  We  choose  Av  so 
that  Av<0  .  The  relationship  between  these  parameters  is 

*  *  *  * 
v  v  +Av  r 

—  *  — T7 —  ■  —  =  cos  <J> 
v  v  +Av  r 


(3.4) 


(3.5) 


-  10  - 


Moreover,  (3.4)  and  (3.5)  imply  the  inequalities 

u£v<v+AvSw  , 

and  (3.6) 

* 

Av  <  6  cos  <f> 

Spheres  centered  within  this  layer  generate  on  M  disks  with  centers 
between  two  parallel  lines,  whose  distances  from  0  are  v-6  and  v+Av+8  . 
Furthermore,  centers  of  disks  which  fall  between  the  parallel  lines  of 
distances  v+8  and  v+Av-6  follow  a  homogeneous  Poisson  process.  The 
intensity  of  the  homogeneous  Poisson  process  for  each  such  layer  is  pAv  , 
where 


■»  U 

“  - 1  ?  • 


(3.7) 


Consider  a  sphere  in  the  prism  P  ,  centered  on  the  plane  1/  in 
distance  d  from  M  .  This  sphere  generates  on  M  a  disk  whose  center 
falls  on  a  line  of  distance  z  from  0  (see  Fig.  2).  Moreover, 

Iz-vl  "  d@/b 

Accordingly,  if  x  is  the  radius  of  the  sphere,  the  radius  of  the  disk 


(3.8) 


y  -  f(x2  -  (z-v)2b2/S2) 


Let  G(y|z,v)  denote  the  conditional  c.d.f.  of  the  random  radius  of  a 
disk,  Y  ,  given  that  its  center  is  on  a  line  having  distance  z  from 
0  and  the  center  of  its  generating  sphere  is  on  a  plane  intersecting 
M  at  a  distance  v  from  0  .  According  to  (3.9),  this  conditional 
c.d.f.  is  given  by 


G(y|z,v)  -  F^(y2  +  (z-v)2b2/62)1/2N| 


It  is  interesting  to  note  that 


I  ■ 


(3.9) 


(3.10) 


G(0 1 z,v)  -  F(  £  I z~v j \ 

'  6  i 

Furthermore,  G(y|z,v)  *  l  for  all  y  >  b(l  -  (z-v)2/g2) 


(3.11) 
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4.  Moments  of  The  Visibility  Measure  on 


In  the  present  section  we  define  a  measure  of  visibility  on  an  interval 
7  on  C  and  its  moments.  This  definition  and  the  derivation  of  the  visi¬ 
bility  moments  on  7  follow  the  general  methodology  of  [  4  } ,  which  was 
developed  for  the  two-dimensional  case.  Since  the  curve  C  and  the  bound¬ 
aries  U  and  W  are  parallel  lines,  it  is  more  convenient  to  formulate  the 
various  formulae  in  terms  of  cartesian  coordinates,  rather  than  using  the 
general  formulae  of  [4],  which  are  based  on  polar  coordinates.  In  Fig.  4, 
we  present  the  geometry  of  the  shadowing  process  on  M  .  The  source  of  light, 
0  ,  is  the  origin  of  the  coordinate  system.  The  x-axis  is  parallel  to  the 
line  C  .  The  centers  of  the  disks  are  distributed  in  a  ■  *ip,  S  ,  between 
two  parallel  lines,  whose  Z-coordinates  are  u-S  and  \  ,  respectively. 

The  interval  of  interest,  C  ,  is  bounded  by  the  end  poin  (s',r)  and 

(s'^r)  .  A  point  P  ■(s.r)  in  C  is  called  visible,  i  ke  line  segment 

s 

OP  does  not  intersect  any  shadowing  disk.  The  measure  o  jibility ,  V{C}  , 


is  the  total  portion  of  segments  of  C  which  are  visible.  Let 


,  if  Px  ■  (x,r)  is  visible 


,  otherwise 


(4.1) 


/  I(x)dx 


(4.2) 


If  p(s)  is  the  probability  that  the  point  Pg  =  (s,r)  is  visible,  then 


r  -  s' 
:iV{C}  /-  / 


p(s)ds 


(4.3) 


Similarly,  if  s^,...,sn  are  the  x-coordinates  of  n  points  in  C  ,  and 
p(s^, ...,sn)  is  the  probability  that  these  n  points  are  simultaneously 
visible,  then,  as  shown  in  [4],  the  n-th  moment  of  V{C}  is  given  by 

E-V  {C}/-  *  n  if  ....  /  p(s,,...,s  )  ds-...ds 
.  1  n  l  n 

A 


where  A  =  {s'^s.  s  Ss''} 
l  n 


(4.4) 
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The  visibility  probabilities  p(s^, . . . .s^)  ,  for  all  nil  ,  are  given,  as 
in  [  4 ]  ,  by  employing  the  functions  defined  below. 

Let  AK  (s,t)  be  the  expected  number  of  disks  falling  within  the  inter¬ 
section  of  the  triangle  0ps_cps  anc*  the  strip  S  ,  which  do  not  inter¬ 
sect  the  line  OP  .  (We  include  also  disks  with  zero  radius).  Similarly, 
s 

we  define  the  function  AK+(s,t)  ,  which  corresponds  to  the  triangle  0psps+t 
(see  Fig.  4).  These  functions  will  be  explicitly  developed  for  the  standard 
(B=0)  and  non-standard  (S>0)  cases  in  the  following  sections.  It  will  be 
shown  that  K  (s,t)  ■  K+(s,t)  *  K(s,t)  .  Let  C  be  a  subset  of  the  strip 
S  that  contains  the  centers  of  all  disks  which  may  cast  shadows  ou  C  (see 
Fig.  3).  This  set  is  obtained  by  the  intersection  of  S  with  the  triangle 
op6'p9t«  »  where 


0'  <  s'  -  b  ((s')2  +  r2 


and 


1/2 


/(u-e) 


o 


2  2\1/2 

9”  >  s"  +  b  f(s’*r  +  r  /  /(u-3) 


(< 


Let  v{ C)  be  the  expected  number  of  disks  falling  in  C  (including  ones 
having  a  zero  radius) .  Explicit  formulae  for  the  standard  and  non-standard 
cases  will  be  given  later.  Accordingly,  the  expected  number  of  disks  that 


may  cast  shadow  on  a  point  Pg  in  C  is  v{C)  -  A(K(s,s-9')  +  K(s,9,,-s)) 


From  the  Poisson  field  assumptions  we  obtain  that  the  probability  that  P 


is  visible  is 


p(s)  *  exp [ v { C}  -  a(k(s,s-9')  +  K  (s,9,,-s)'j 


\ 


0 


For  any  n  points  P  ,...,P  in  C  ,  such  that  s'<s_<...<s  <s' 


the 


1  n 

probability  of  their  simultaneous  visibility  is 


r 


p(s1,...,sn)  =  exps-[v{C,  -  A  |K(s^,s^-9 ' )  + 


0 


K(sn’9"-sn))  '  A  (K(si’VSi) 


+  K(Si+i,si+i-s.)!  ] 


where 


,i  / _ -iri 


if  (can 


'S  . 


+  tan 


-1 ,  Si+1 


r  tan 


( 


The  point  P~  is  chosen  on  the  bisector  of  the  angle  between  OP 

Si  Si 

and  OP  in  order  to  assure  that  a  disk  whose  center  is  in  the 

Si+1 

sector  between  OP  and  0P~  ,  which  does  not  intersect  OP  ,  will 

Si  si  si 

not  intersect  OP  ,  and  vice  versa. 
si+l 

The  probability  that  the  whole  interval  C  is  visible  is 

PL  -  exp  -tv*{c}  -  X  KCs'.e'-e')  +  K(s,,,0',-s")  ] 

£ 

where  v  {C}  is  the  expected  number  of  disks  falling  in  S  between 

OPg ,  and  0PQ , ,  ,  minus  the  expected  number  of  disks  falling  in  S  , 

between  OP  ,  and  OP  , ,  having  zero  radius . 
s  s 

Finally,  as  shown  in  [  3]  and  [  4 ]  ,  the  moments  of  V{5}  can  be 
determined  recursively,  in  the  following  manner.  Let 

’Pq(s)  »  exp{XK(s,s-6') }  ,  s'<s^s" 

and,  for  every  j>l  ,  define  recursively 


'Pj  <s)  -  /  'J'j^Cy)  expiAtKCy.^Cs))  +  K(s ,y2  (s) )  ]  }  dy 


where 


and 


y1(s)  =  r  tan 


-  y 


y2(s)  -  s-y-y^(s) 

The  n-th  moment  of  V{C}  is  given  then  by  the  formula 


y  *  n!  exp{-v{C}}  J  y  .(s)  exp{AK(s,9 ' ’-s) }  ds 

H  ■  11“  J. 

a  * 


(4.10) 

(4.11) 

(4.12) 

(4.13) 

(4.14) 

(4.15) 


In  the  following  section  we  develop  more  explicit  formulae  for  these  functions. 


5.  More  Explicit  Formulae  for  K(s.t)  and  v{C} 


5.1  The  Standard  Case  (3*0) 

In  the  standard  case,  the  expected  number  of  disks  in  C  is  the  intensity, 
2Ab  ,  multiplied  by  the  area  of  C  ,  i.e., 

a  *  «_o l  2  2 

v{C}  =  Ab  ~~  (w  -u  )  , 

and  (5.1) 

v*{C}  »  v{C}  -  ||  (w2-u2)(s"-s•)  . 

The  geometry  of  the  standard  case,  (5.1)  as  illustrated  in  Fig.  3,  shows  that 
the  location  of  centers  of  disks  of  radius  y  ,  which  are  between  the  rays 
OP  and  OP  ,  and  which  do  not  intersect  the  ray  OP  ,  is  the  trapezoid 

S**u  S  S 

EyFyCyDy  .  Similarly,  the  trapezoid  Ey  Fy  Dy  contains  the  centers  of 

all  disks  of  radius  y  ,  which  are  between  OP  and  OP  , .  ,  which  do  not 
_  1  *  s  s+t 

intersect  OP  .  Notice  that  the  area  of  E  F  C  D  is  the  same  as  that  of 
s  y  y  y  y 

EyFyCyDy  .  Finally,  since  K_(s,t)  and  K+(s,t)  are  the  average  areas  of 
these  trapezoids,  with  respect  to  the  distribution  G(y)  ,  multiplied  by 
2b  ,  we  obtain  that 

K_(s,t)  ■  K+(s,t)  2  K(s,t)  ,  for  all  s'^sSs”  (5.2) 

and  all  0  <  t  . 

We  develop  now  the  formula  for  K(s,t)  .  Notice  first  that  the  area  of  the 

trapezoid  EyFyCyDy  can  be  computed  as  the  area  of  the  triangle  ByCyDy 

minus  the  area  of  the  triangle  ByEyFy  .  Let  A(s,t,v,y)  be  the  area  of 

such  a  triangle,  whose  vertices  are  B  C'D"  ,  where  C"  and  D"  are 

y  y  y  y  y 

on  a  line  parallel  to  the  x-axis,  at  distance  v  from  the  origin.  The 

triangle  B  C^D'1  is  similar  to  the  triangle  OP  .P  ,  since  B  D*  I  I  OP 
_ 2LZ  y  w  s-t  s  y  y  11  s 

Hence,  C''D''+d  *  —  t  ,  and  the  area  of  the  triangle  0C,'A''  is 

2  y  y  r  y 

•  The  area  of  the  parallelogram  ByA'Dy'A''  is  d(v-h)  .  The  area  of 

the  triangle  OByA'  is  dh/2  .  Moreover,  h  =  |  d  ,  and  since  the  triangles 

2  2  1/2 

AP  ,P  and  OP  P  are  similar,  d  *  y(s  +r  )  / r  .  From  all  these  geomet- 

S"*U  S  €1  S 

rical  considerations  we  obtain  that. 
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,  .  2.  2,1/2  ,2  s2+r2 

(vt/(s  +r  )  -y)  2tr 


if  y  ^  vt/(s2+r2)1/2 


A(s,t,v,y)  H 


(5.3) 


,  otherwise 


Finally,  let 


s2+r2  mifn(b,vt/(s2+r2)1/2)  /  2  2  1/2  N2 

K(s,t,v,r)  =  b  — —  /  ^y  -  vt/(sZ+r  )X/  j  dG(y)  ,  (5.4) 

A  A 

then  K(s , t)  *  K(s,t,w,r)  -  K(s,t,u,r) 

In  Section  6  we  further  develop  this  formula,  for  the  uniform  distribution 
F(x)  .  In  this  special  case  G(y)  is  given  by  formula  (3.3). 

5.2  The  Non-Standard  Case  (g>0) 

The  expected  number  of  disks  centered  in  C  is,  in  the  non-standard 
case,  the  same  as  in  the  standard  case.  Indeed, 


Ab  e"-e'  7  a  V*  . 
v{C}  ■  —  •  -  J  dv  J  z  dz 

P  r  u  v-g 

-  (e"-0')(w2-u2) 


(5.5) 


In  a  similar  fashion  we  obtain  that 


,,  , ,  ,  w  v+g 

Ab  s'-s'  r  .  r 


'{C}  *  v{C)  -  —  •  ^ —  /  dv  /  z  F ( 


u  v-g 


W  X 

v{C)  -  —  (s  *  *  — s 1 )  j  dv  f  (u+vg)  F(b|u|)  du 

u  -1 


In  the  special  case  of  uniform  distribution 


F(x)  ■  r-  min(x,b)  ,  OSx  ,  one  gets. 


v*{C}  -  v(C}  -  £  (s”-s,)(w2-u2) 


The  K-functions  are  obtained  in  the  following  manner.  As  seen  in  Fig.  4, 

consider  a  disk  with  a  center  on  a  parallel  line,  Z  ,  intersecting  the 

axis  at  z  ,  u<-B£zSw-H3  ,  and  having  distance  x  from  the  ^intersection  of  Z 

with  0P„  .  The  maximal  radius  that  such  a  disk  can  have,  without  inter 

s _ 

secting  OP  ,  is 
s 


y(x) 


xr/(s2+r2) 


1/2 


.  z 
Osxs -  t 
r 


(5.6) 


Since  this  structure  is  the  same  for  disks  on  the  right  or  the  left  sides 
of  OP  ,  K  (s,t)  ■  K  (s,t)  =  K(s , t)  ,  where 

S  "m  » 


,  w  v+e  Zt/r  ,  -  1/2 

K(s,t)  *  t  /  dv  /  dz  J  G  (xr/(s  +r  )  I  z,v)  dx  (5.7) 

p  u  v-0  0 

Substituting  (3.10)  in  (5.6)  and  making  simple  change  of  variables,  one 
obtains 

2  2  ^2 

.  ,  ,l/2w  1  fvffjz)t/(s%r  )  ,  ,  ,  1/2 

K(s,t)  «  -  (s^+r  )  /  dv  /  dz  /  F((y  ■+*  O  )  dy 


.  (5.7) 


Fig.  4.  The  Geometry  On  M  :  Non-Standard  Case. 
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6.  Further  Development  for  the  Uniform  Case  and  Numerical  Illustrations 


6.1  The  Standard  Case 

In  the  present  section  we  provide  more  explicit  expression  for  the  function 

A 

K(s,t,v,r)  ,  given  in  (5.3),  when  the  distribution  of  the  radii  of  spheres, 

F(x)  ,  is  uniform  on  (0,b)  ,  i.e.  F(x)  -  ^  min(x,b)  ,  for  all  0<x<®  .  The 

corresponding  formula  of  G(y)  ,  for  the  standard  case  if  given  in  (3.3).  i/o 

2  2\ 

We  distinguish  between  two  cases:  Case  I,  the  case  of  b  2  vt//s  +r  J 

2  2  1/2 

and  Case  II,  the  case  of  b  >  vt/(s  +r  ) 


Case  I 


According  to  (5.3)  , 


2  2  j 

K(s,t,v,r)  -  b  E{Y2}  -  2 


2. _2\l/2 


2  2 
s  +r 


(6.1) 


Formulae  for  the  first  two  moments  of  Y  in  the  standard  case  are  developed 
in  Appendix  A.l.  According  to  these  formulae 

2  2 

K(s,t,v,r)  -  .291667b3  -  .745922b2  v(s2+rV/2/r 


+  bt  v  /r 


(6.2) 


Case  II 


In  Case  II, 


2  2  ? 
s  +r  ,  ? 


K(s,t,v,r)  *  — — -  b  /  (y-£)  d  G(y) 


where 


Equivalently, 


C  -  vt/  ^s2+r2j 


«  2  Ly  2  2  £ 

K(s,t,v,r)  *  bv  t/r  -  2r~*  s  +r  *  J  y  dG(v) 

0 

+  s+L-  b  J  y  2  dG(y) 


(6.3) 


(6.4) 
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In  Appendix  II  we  evaluate  the  incomplete  moments  ECY1;?}  ■  /y^dG(y)  ,  i»l,2  . 

0 


According  to  (6.4)  and  formulae  (A. 12)  and  (A. 19)  we  obtain 


2  r 

K(s,t,v,r)  -  bv2t/r  -  ^^(s^r2)1^2  •  ^  ~  tan"1 


imi 


^  1/2 


,2  I 


-  ios<H V  -R>  ii  -  <£>2')1/2 


\7'b''6  /  2"'b'  YVC 


+  (£)3  log  I  1  +  <1  -  <j)2)1/2J  b3  <ji  - 

i  (l  -  (£>2)  1/2  (l  +  i<£>2)  -  i(f)‘  log  ('  <{b)/(l 

a  -  tff'h  )(. 


where  A(z)  -  .0933z2  +  .57z^  -  .3633z6 


(6.5) 


6.2  The  Non-Standard  Case  (8>0) 

The  function  K(s,t)  in  the  non-standard  case  is  evaluated,  for 
the  uniform  distribution  of  sphere  radii,  but  substituting  F(x)  *  ^  min(x.b)  , 
0Sx<«  ,  in  formula  (5.7).  Since 

((£)2+z2)1/2  ,  if  y  s  b(l-z2)1/2 


r  ,  2^.  2  2.1/2 

F  (y  +b  z  ) 


,  if  y  >  b(l-z2)1/2 


we  distinguish  between  two  cases: 

Case  I  (v+z8)t/b(s2+r2)1/2  s  (l-z2)1/2  . 

2  2  1/2 

Let  y(t)  -  t/(s  +r  )  '  and  g(z,v,t>  »  y(t)(v+z0)/b 


Then,  In  Case  I, 
b|(z,v,b) 

0 


|(2,v,b)  ^  ^2^2^!/^  ,  b  |g(8fVtt)  (g2(ZfV>t)  +  z2)  i/2 


+  Z2  log (g (z , v, t)  +  g2 (z , V, t)  +  z2)1/2 


)  - 


z2  log  Iz 1 1 


(6.7) 
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Case  II  (v  +  z8)t/b(s2+r2)1/2  >  (l-z2)1/2 


In  Case  II, 


b*(z,v,t)  2  2  2.1/2,,  b  f2,  .....  2,1/2, 

J  F((y  4b  z)  '  )dy  -7|Z  log(l+(l-z  )  ) 

n  z  L 


-  z  loglzl  +  2g(z,v,t)  -  (1 


-*2>1/2~!  . 


Thus,  from  (5.7),  (6.7)  and  (6.8)  , 


(6.8) 


,2  ,  2^  2,1/2  w  1  ,  ,/0 

K(s, t)  -  2 - - - p -  J  dv  /  dz  jl{g(z,v,t)  S  (l-z  )  } 

u  -l  L 


u  -1 


^g(z,v 


.t)(g  (z,v,t) 


2Y/Z  +  z2log(g(z,v,t)  +  ;g2(z,v,t)  +  z 


+ 


-z2log|zA  +  I{g(z,v,t)  >  (l-z2) 1/2 }  (z2log(l+(l-z2)1/f2  j 


-z2log|z|  +  2Y(z,v,t)  -  (1-Z2) 1/2 y  j  ,  (6.9) 

where  I(A>  designates  the  indicator  function  of  the  set  A  .  The 
double  Integral  In  (6.9)  was  evaluated  numerically,  applying  an  iterative 
procedure  with  a  20-point  Gaussian  quadrature  (see  Abramovitz  and  Stegun 
[1;PP.  916]). 


6.3  Numerical  Illustrations 

In  the  present  section  we  provide  some  numerical  computations ,  which 
illustrate  the  above  theory.  The  following  parameters  are  held  fixed 
throughout  the  various  cases. 

1.  The  distance  of  U*  from  the  origin  is  u*  ■  .4  . 

2.  The  distance  of  W*  from  the  origin  is  w*  ■  .6  . 

3.  The  distance  of  C*  from  the  origin  is  r*  •  .1  . 

The  radii  of  the  spheres  are  i.i.d.  random  variables  having  a 

uniform  distribution  on  [0,b]  ,  with  b  ■  .3  . 

The  line  segment  ^  is  of  2  units  length,  and  is  located  between 

s’  ■  -1.  and  s”  ■  1.  The  varying  parameters  are:  the  Intensity  of  the 

Poisson  field,  X  m  2,6,10  and  the  inclination  angle  of  the  plane  M  , 

<P  m  0  ,  ir/8  ,  tt/4  and  ir/3  .  In  Table  6.1  we  present,  for  each  one  of 

these  twelve  cases,  the  first  ten  standardized  moments  of  V{C>  and 

the  moments  of  the  approximating  mixed-beta  distribution.  The  standardized 

moments  of  V are  u  ■  y  /(s,,-s,)n  ,  n  *  1,2,...  ,  where  p  have 

n  n  n 


been  determined  according  to  the  recursive  formula  (4.19).  The  mixed-beta 
distribution  is  a  distribution  function  on  [0,1]  ,  having  a  c.d.f. 


y0 +  a  -  po  -  P1>  T?(6r  l  "c“1  <1-)B'1  *■  •  0M<1 

F*(x)  (6 

t  •  ls*  ' 

This  is  a  beta-distribution  on  (0,1)  ,  mixed  with  a  two-point  distribution, 
with  jumps  of  height  Pq  at  x*0  and  at  x»l  . 

P0  -  P{V{C}  -  0}  and  P^^  -  P{V{C)  -  s,,-s*} 

P^  is  given  by  formula  (4.14).  In  [3]  and  [4]  we  have  shown  how  the 

parameters  P_,  a  and  8  can  be  determined,  by  equating  the  first  three 
”  _  * 

standardized  moments  of  V{C}  ,  u  ,  n*l,2,3  ,  to  the  first  three  moments 
*  n 

of  F  (x)  .  In  Table  6.2  we  present  the  parameters  of  the  approximating 

* 

mixed-beta  distribution,  F  (x)  ,  for  the  twelve  cases  under  consideration. 
The  method  of  computing  the  convolutions  in  the  iterative  formulae  of  the 
moments  was  described  in  [4].  To  apply  those  formulae  to  the  present  cases 

ft  £ 

we  substitute,  r-r  /cos  $  ,  w«w  /cos  $  ,  y»y  /cos  <j>  and  fj-b  tan  <J>  ,  in 

the  proper  places.  We  see  in  Table  6.1  and  Table  6.2  that  for  each  intensity 

X  ,  the  total  visibility  probabilities,  P^  ,  as  well  as  the  moments  of  V{C> 
decrease  as  the  Inclination  angle  <t>  increases. 

Table  6.1  shows  again  very  good  agreement  between  the  first  ten  standardized 
moments  of  V{7?}  and  those  of  the  approximating  beta-mixture  distribution. 
Thus,  one  can  apply  the  beta-mixture  distribution  to  approximate  the  c.d.f. 
of  V{7}/(s' ’-s')  .  Thus, 

PQ  +  (l-P0-P1)Iv(ot,8)  ,  0sv<l 

P{V{C}  <;  v(s,,-s')}  - 

1  ,  v>l 

where  I  (a, 8)  is  the  Incomplete  beta  function  ratio. 


A 

D 

2 

3 

n 

5 

6 

B 

8 

9 

10 

QO 

2 

0 

.958 

.932 

.914 

.902 

.894 

.889 

.886 

.886 

.887 

.887 

.841 

.958 

.932 

.914 

.902 

.893 

.886 

.881 

.876 

.873 

.870 

.841 

n/8 

.955 

.927 

.909 

.896 

.888 

.882 

.879 

.878 

.879 

.882 

.831 

.955 

.927 

.909 

.896 

.886 

.879 

.873 

.868 

.865 

.861 

.831 

ir/4 

.944 

.910 

.887 

.871 

.860 

.853 

.848 

.846 

.84S 

.848 

.791 

.944 

.910 

.887 

.871 

.859 

.850 

.843 

.837 

.832 

.828 

.791 

3tt/8 

.904 

.847 

.809 

.783 

.765 

.752 

.743 

.737 

.733 

.732 

.678 

.904 

.847 

.809 

.783 

.764 

.749 

.783 

.729 

.721 

.715 

.678 

6 

0 

.878 

.809 

.765 

9 

.713 

.698 

.688 

.680 

.675 

.674 

.594 

.878 

.809 

.765 

.712 

.696 

.683 

.672 

.664 

.654 

.594 

ir/8 

.871 

.798 

.751 

.720 

.697 

.681 

.670 

.662 

.657 

.655 

.573 

.871 

.798 

.751 

.719 

.696 

.679 

.665 

.654 

.646 

.639 

.573 

tt/4 

.842 

.754 

.699 

.661 

.635 

.616 

.602 

.593 

.586 

.582 

.494 

.842 

.754 

.699 

.661 

.634 

.614 

.598 

.586 

.576 

.567 

.494 

3tt/8 

.739 

.609 

.532 

.482 

.448 

.423 

.406 

.392 

.382 

.375 

.284 

.739 

.609 

.532 

.482 

.447 

.422 

.402 

.387 

.375 

.365 

.284 

10 

0 

.806 

.703 

.640 

.599 

.570 

.549 

.534 

.523 

.515 

.510 

.420 

.806 

.703 

.640 

.599 

.569 

.547 

.530 

.517 

.506 

.497 

.420 

tt/8 

.794 

.687 

.622 

.579 

.549 

.527 

.511 

.500 

.492 

.486 

.396 

.794 

.687 

.622 

.578 

.548 

.525 

.508 

.494 

.483 

.474 

.396 

tr/4 

.750 

.626 

.551 

.503 

.470 

.446 

.429 

.416 

.407 

.309 

.750 

.626 

.551 

.503 

.469 

.444 

.425 

.411 

.309 

3tt/8 

.603 

.440 

.352 

.299 

.264 

.240 

.223 

.210 

.123 

.603 

.440 

.352 

.299 

.263 

.238 

.220 

.206 

.195 

.987 

.123 

*)_  _ 

The  upper  value  in  each  cell  is  the  standardized  moment  of  V{C}  .  The  lower 

value  is  that  of  the  beta-mixture  distribution .  y*»P,  . 

oo  i 
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Appendix  I.  The  Moments  of  Y  in  the  Standard  Uniform  Case 
b 

Writing  E{Y}  *  J  (1-G(y))dy  and  substituting  (3.3)  we  obtain 
0 

E{Y}  =  |  J  (l-(J)2)1/2dy  +  |  /  (£)2log(£)dy 

-  \  J  (^)2log(l  +  /l-(£)2)  dy  •  (A.l) 

Making  the  changes  of  variable,  u=y/b  ,  we  obtain 

E{Y}  =  /  (l-u2)1/2du  +  /  u2  log  u  du 

Z  0  0 

-  /  u2log(l  +  /l-y2)du)  .  (A. 2) 

0 

2 

Furthermore,  making  the  transformation  w*u  ,  we  obtain 

J  (l-u2)1/2du  =  j  j  w"1/2(l-w)1/2dw 
0  1  0 

=  |fi(y,|)  =  */4  ,  (A.  3) 

where  B(a,b)  designates  the  beta-function  at  a  and  b  .  In  addition  the 
other  integrals  of  (A. 2)  are  (see  H.B.  Dwight  [2;pp.  141]) 

1  ? 

J  u  log  u  du  =  -1/9  , 

0 

/  u2log(l  +  /l-u2)du  -  .1505865  ... 

0 

Substituting  these  results  in  (A. 2)  we  obtain  the  first  moment 
E{Y)  *  b{|  -  Ta  -  |(. 1505865)}  =  .261850  b 


(A. 4) 


The  second  moment  of  Y  is  obtained  in  a  similar  fashion  as 


E{Y2} 


h* 

K-' 

L' 

It-' 

K' 

K' 


B 


2  J  y [l-G(y) ]dy 
0 

b  J  (?)a-(J)2)1/2dy  +  b  /  (£)3log(£)dy 


-  b  /  (£)3log(l  +  (1  -  (^)2)1/2)dy 

-  .166667 


Appendix  II.  The  Incomplete  Moments  of  Y  In  The  Standard  Uniform  Case 
The  first  two  incomplete  moments  of  Y  are 

.  5  i 

E{Y  }  =  /  y  dG(y)  ,  i»l,2  ,  where  £<b 

0 

The  incomplete  moment  of  the  first  order  is 


E{Y;?}  =  J  [1-G(y)]dy  -  £[1-GU)] 
0 


(A. 5) 


(A. 6) 


(A.  7) 


Substituting  (3.3)  for  G(y)  we  obtain  that 


/  (1-G(y))dy 
0 


|  J  (l-(J)2)1/2dy  +  j  f  (£)2log(£)dy 
0  0 


-  j  J  (^)2log(l  +  /l-(£)2)dy 


Using  a  result  from  Abramowitz  and  Stegun  [l;pp.  946],  we  can  show  that 


1  S  /■,  ,y* 1/2  ,  br  -1,  (c/b)‘  *1/2  .  c  „  ,£*2*1/2, 

j  J  (1  -  (p  )  djr—  ^[tan  (  w  ■■■  2>  +  ^  (!  -  (^)  )  ]  .  (A. 9) 


(A.  8) 


Furthermore 


|  |  (£>2log(f)dy  -  |(f)3Uog(f)-i]  . 

2  1/2 

Finally,  by  making  the  transformation  w=(l-(y/b)  )  ,  we  obtain 

2 

J(£)2log(l  +  (l-(y/b)2)1/2)dy  -  \  (?Jb)  w1/2log(l  +  (l-w)1/2)dw 
0  ^  0 

-  |[.0933(^)2  +  .57(£)4  -  .3633(£)6]  . 


(A.  10) 


(A. 11) 


The  approximation  on  the  RHS  of  (A. 11)  was  obtained  by  approximating 
1/2  1/2 

w  log(l  +  (1-w)  )  on  the  interval  (0,1)  by  the  2nd  degree  polynomial 

.0933  +  1.14w  -  1.09w2  .  From  (A. 7)  -  (A. 11)  and  (3.3)  one  obtains, 

E{Y;?}  -  b{^tan~1(-(g/b)22-  )1/2  -  ±(£)3(log(£)  +  ±) 

4  l-(5/b)2  3  b  b  6 

-  -  £<£><1  -  (£>2)1/2  +  Y(£)3log(l  +  (1 

where 

A(z)  =  .0933S2  +  .57z4  -  .3633z6  .  (A. 13) 

The  incomplete  moment  of  the  second  order  is 

E{Y2; 5}  -2  /  y[l-G(y)]dy  -  e2(l-G(0)  •  (A.14) 

0 


(A. 12) 


-  (£)2)1/2)> 


From  (3.3)  we  obtain 

^  £  £ 

2  /  y[l-G(y)]dy  -  /  y(l-(^-)2)1^2dy  -  b  /  (£)  3log (£) dy 
0  0  0 

-  /  (f)3log(l  +  (l-(£)2)1/2)dy 


(A. 15) 


-  29  - 


Following  Che  integration  methods  applied  above,  we  obtain 


and 


^-[l-(l-(f)2)3/2] 


(A.  16) 


S  2 

b  j!)  ^3l°8^dy  "  V  *  ^b^log^  “  ^  *  (A. 17) 

2  1/2 

Furthermore,  by  making  the  transformation  w  *  1  +  (l-(y/b)  )  ,  we  obtain 

b  /  (£)3log(l  +  (l-(£)2)1/2)dy  - 
0  b  b 

2 

b 2  /  2  ,  .A  (w-l)-(w-l)3]log  w  dw  =  i  -  -i-(l-(^-)2)1/,2(l  + 

l+(l-(?/b)2)1/2  6  6  b 

i(b)2)  +  H)4(loga-Kl-(f)2>1/2)  -  (A.  18) 

Finally,  substituting  (A. 15)  -  (A. 18)  in  (A. 14)  we  arrive  at  the  formula 

E{Y2;C}  -  b2{|  -  |(1-(£)2)1/2(1  +  |(^)2) 

-  \  log((?/b)/(l  +  (1-(|)2)1/2)}  .  (A. 19) 
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